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ABSTRACT 

In the present paper, an overview of the aerodynamic aircraft-design methods and their 
recent applications in Japan is presented. One of the design codes which was developed at the 
National Aerospace Laboratory (NAL) and widely used now is mainly discussed, and hence, 
most of the application examples are the results of the collaborative works between heavy 
industries and National Aerospace Laboratory. Wide variety of applications in transonic to 
supersonic flow regimes are presented. Although design of aircraft elements for external flows 
are the main focus, some of the internal flow applications arc also presented. Recent applications 
of the design code using the Navier-Stokes and Euler equations in the analysis mode include the 
design of HOPE(space vehicle) and USB(upper surface blowing) aircraft configurations. 


INTRODUCTION 

With the advent of supercomputers having fast processors and large memories, 
CFD(computational Fluid Dynamics) is progressing at incredible speed. Three-dimensional 
Navier-Stokes simulations, which were very rare ten years ago even for relatively simple body 
configurations are now common at any conference on fluid dynamics! 1,2, 3]. Flow field 
simulations over complex body configuration are not difficult task once the geometry data is 
given. We can learn a lot of flow physics from the simulated results that may be helpful for re- 
designing the body configuration. Although such simulated results give us a lot of information 
about the flow field, they would not tell us how to modify the body configuration for the better 
design. One way to do it may be a trial-and-error type approach where conducting a large 
number of simulations is necessary, which is still not feasible even with advanced 
supercomputers. So-called design programs for determining the optimum geometry may be as 
useful as analysis programs simulating the given flow fields. 

There has been a strong effort to develop both airfoil and wing design methods for many 
years. Unfortunately, the progress is not as remarkable as analysis methods.This is true in Japan 
as well as in the United States. CFD technology has been remarkably improved last several years, 
but on the other hand, no much progress was made for the design methods and code 
development. Only one remarkable progress in Japan was the design method developed by 
Takanashi at National Aerospace Laboratory in 1984. His method is "iterative correction method" 
I based on the perturbation equations of potential flows. In this method, the geometry correction is 
L made iteratively to reduce the difference between the target pressure distributions and the 
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r pressure distributions that is obtained by an analysis code. One of the advantages of this code is 
that any analysis code can be incorporated because analysis code is sort of black box for the 
geometry correction method. Analysis codes are not necessarily restricted to the potential codes. 
Even the Navier-Stokes codes can be used although the convergence is not guaranteed. Because 
of the flexibility and robustness of the code, it has been used for wide variety of applications. 
Now, most of the aircraft industries in Japan use this computer code and applied it to the practical 
problems. 

In the present paper, Takanashi's design method and its applications are presented. Since 
this is a paper giving an overview of the Japanese activity, only the conceptual explanation is 
given about the method itself, and the focus is mainly laid on the demonstration of the 
applications to a wide variety of the flow fields. 


DESIGN METHODOLOGY 


Background 

There are several approaches for the design problems. One way may be the numerical 
optimization using an analysis code. Wing design method was proposed by Hicks in 1976[4], 
and the research has been extensively conducted since then. In this optimization technique, a 
wing section with, for instance, minimal total drag under some constraints such as a specified lift 
and maximum thickness is sought by using the analysis code and the optimization code 
iteratively. Recently, Jameson[5] proposed an efficient method using a control theory. There 
exists so-called "inverse method" of wing design in which wing geometry is determined to realize 
the specified pressure distributions. This type of approach was used for wing design by 
Henne[6] and for wing-fuselage design by Shankar[7] for example. The approach used by 
Takanashi may be different from either of the approaches above. This is an iterative residual 
correction method similar to the works by Barger and Brooks[8], Davis[9], and McFadden[10] 
for the two-dimensional problems. The advantage of this approach is that only minimum effort in 
developing the geometry correction code is needed to decrease the pressure residual, while an 
analysis code is retained in its original form. In the next section, the formulation is briefly 
described. 


Formulation of Inverse Problem and Iterative Procedure 

Only a concept of the design method that was developed by Takanashi in 1984 is briefly 
described. More details can be found in his original and the following papers[l 1,12,13]. 

First, inverse problem is defined. Here the nonlinear full potential equations are taken as 
basic equations, and in the formulation process, small perturbations are assumed. Thus, the 
applicability is restricted to the flow field without shock waves or with weak shock waves. After 
some manipulations, integral equations that relate the geometry change and the surface pressure 
change are formulated. Iterative design procedure is formulated using the integral equations 
obtained above. Body (wing, wing-body complete aircraft etc.) surface is paneled into segments 
and the integral equations are discretized and numerically solved to find the necessary amount of 
geometry modification once the difference of required and calculated pressure difference is 
defined. Since we have the target pressure distributions which is required, we can define the 
difference using some analysis code. 

The iteration process can be defined as follows. First, we assume initial body geometry, 
then calculate the surface pressure distributions using some analysis codes. Since we know the 
|_ require pressure distributions, we can calculate the difference between the required an calculated 
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r tip to satisfy the off-design requirements is shown in Fig. 9. The wind tunnel experiment was 
conducted to check the aerodynamic performance of the designed wing. The measured Cp 
distributions are presented in Fig. 10. Reasonable agreement is observed between the target and 
the measured pressure distributions. Figure 1 1 shows the comparison of the pressure contours on 
the upper surface of the wing. Here again, good agreement is obtained between the computed and 
measured contour plots even though small discrepancy is observed near the root and tip. 

At the time of this design code development, there was a collaboration between JADC 
representing Japanese industries and Boeing company to develop a new transonic aircraft. The 
project was called 7J7 in the United States, and YXX in Japan. Although this project was retailed 
because of the market change, there left is a lot of technology accumulations for die Research and 
Development Under this project, many wing configurations were designed by Mitsubishi Heavy 
Industries again using Takanashi's code. Some of the designed wings were used for the 
simulations using the Reynolds-averaged Navier-Stokes equations! 17, 18] and the computed 
results were compared with the corresponding experiments! 19] to confirm the aerodynamic 
performance of the designed wings. These examples will be shown at the conference. 


Airfoil Design Using Navier-Stokes Equations 

As has been mentioned above, the analysis code is sort of a "black box" and it can use any 
analytical method even though the convergence is not necessarily guaranteed. Hirose et al. 
coupled Takanashi's design code with two-dimensional Reynolds-averaged Navier-Stokes 
code[20]. With specifying the same pressure distributions at each spanwise station for large 
aspect ratio wing, the three-dimensional design code was incorporated with the two-dimensional 
Navier-Stokes code for the design of two-dimensional airfoil. One of the application examples is 
shown here. Shockless supercritical pressure distributions at Cl= 0.6 was specified as a target 
and the initial geometry was set up to have strong shock wave. The freestream Mach number is 
0.75 at the Reynolds number 13 million. The initial, target and computed Cp distributions along 
with the initial and final airfoil geometries are plotted in Fig. 12. The target Cp distributions are 
almost realized in ten iterations. 


Two Dimensional Transonic Cascades 

Takanashi reformulated his original design code and developed a two dimensional cascade 
design program in 1986. The analysis code in this case is a Euler code using explicit time 
integration. Even after 10 iterations, fully converged solution was not obtained. However, the 
pressure is becoming closer and closer to the target pressure on every iteration stages. The 
solution after 10 iterations is presented in Fig. 13 along with the cascade geometry. Takanashi 
insisted in his paper! 13] that the convergence would be much improved by optimizing the 
parameters in the design process for cascade flows. 


Additional Applications 

Recently, with the rapid progress of supercomputers, the design code above was 
combined with three-dimensional Navier-Stokes codes and applied to more difficult cases. Both 
Mitsubishi (MHI) and Kawasaki (KHI) Heavy Industries applied it to the design of HOPE( H-II 
Rocket Orbiting Plane). The HOPE is a space vehicle that NASDA (National Space Development 
Agency) is currently developing. Both companies were interested in redesigning the tip fin of the 
configuration. MHI analysed the transonic flow at Mach number 0.9 with 5 degrees angles of 
attack and the Reynolds number 2 million[21]. They found by the Navier-Stokes simulations that 
the flow field surrounded by the fuselage, main wing and tip fin became almost channel flow and 
[_ strong shock wave and associated flow separation occured. The Takanashi's design code was 
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r iteratively used with the three-dimensional Navier-Stokes code and the good improvement was 
obtained after five iterations. They used half a million of grid points for the Navier-Stokes 
analysis and the computer time required for each iteration step was 5 hours for the analysis mode 
and 0.2 hours for the design mode Thus, in total, 26 hours were necessary even with the Fujitsu 
VP400, one of the most advance supercomputers at that time. The initial body configuration is 
shown in Fig. 14 in terms of the computational gnd. The initial and the final chordwise pressure 
distributions on the tip fin are presented in Fig. 15 with the corresponding sectional geometries. 
Although the target pressure is not precisely realized, there is obvious improvement such as 
disappearance of the suction peak. The close-up views of the near-surface streamlines obtained 
from the computed flow fields both for the initial and obtained configurations are presented in 
Fig. 16. Shock wave is weakened and the flow separation on the tip fin surface disappears in the 
final configuration. 

Kawasaki Heavy Industries tried to modify the pressure distributions over the tip fin to 
satisfy the buffet boundary by re-designing the tip fin using Takanashi's design code with the 
Euler code[22]. About 200,000 grid points were used in the analysis mode and total computer 
time for five iterations was about 5 hours. In this example, the freesteam Mach number is 0.9 and 
the angle of attack is 6.5 degrees. The original and designed sectional geometries, and the initial 
and final Cp distributions along with the target Cp are presented in Fig. 17. Remarkable 
improvement is observed although the target Cp distributions are not realized also in this 
example. 

Kawasaki Heavy Industries also applied the design code for the redesign of the USB 
(Upper Surface Blowing) wing configuration of the STOL[23]. The planform of the USB is 
shown in Fig. 18. In this example, Isobar pattern is the target, but the wing section is fixed near 
the nacelle and the tip to avoid resulted very thin wing section to weaken the shock wave. Figure 
19 shows the sectional Cp distributions. The strong shock wave that appeared on the initial 
configuration is weakened and the target Cp distributions are almost realized. 

Another aircraft company named Fuji Heavy Industries developed their own design code 
based on the Takanashi's method. They applied it to the design problem of wing-fuselage 
combi nation [24], The analysis code was full potential code. The target pressure distributions 
were such that realize the isobar pattern on the wing surface and are the same as the initial ones on 
the fuselage. The initial and final Cp distributions and the surface pressure contours are plotted in 
Fig. 20. The computed Cp in the lower surface realizes the target Cp, but still some discrepancy 
exists on the upper surface. However, compared to the initial Cp distributions, improvement is 
obvious. The final configuration is shown in Fig. 21. 


SUMMARY 

An overveiw of the Aerodynamic aircraft-design methods and their recent applications in 
Japan was presented. One of the design codes developed at the National Aerospace Laboratory 
(NAL) is mainly discussed because of its popularity in Japan, and wide variety of applications 
were presented from transonic to supersonic flow regimes. This design method uses inverse 
design code and analysis code iteratively to realize the required pressure disbributions, and thus 
any anaysis code can be used. Some of the examples shown here used Euler and Navier-Stokes 
code as an analysis mode. These application examples indicated the capability and feasibility of 
the design code. The fact that many companies currently use this design code for practical 
problems and obtain successful results proves it. 

This paper is written based on the results that the first author has noticed. There may be 
more activities in Japan that can not be included in the paper. Unfortunately many of the papers in 
the reference list are written in Japanese. However, some of the important papers such as J 
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r Takanashi's original paper are written in English and the authors hope that the list of reference in 
this paper is useful for any researchers for the design problems. 
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Fig. 4 Pressure contour plots on the upper 
surface of the designed wing[l 1]. 

L 



Fig. 3 Designed wing geometry and the 
pressure distributions[ 11]. 



Fig. 5 Designed wing geometry with the plots 
of the pressure distributions and the 
upper surface pressure contour plots 
(potential code and the boundary code 
were used in the analysis mode)[12]. 
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Fig. 6 Target chordwise pressure distributions at each span station[16]. 
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Fig. 8 Initial and final pressu distributions along with the target pressure distributions plotted 
with the wing geometries[ 1 6] . 




Fig. 9 Final wing geometry taking off-design requirements into account! 16]. 



Fig. 10 Experimentally measured and designed pressure distributions! 16], 
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Fig. 1 1 Computed and experimentally measured pressure contour plots on 
the wing[16]. 
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Fig. 14. Computational grid over a HOPE configuration[21]. 
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Fig. 15 Initial and final chordwise pressure distributions over the tip fin of the HOPE and their 
spanwise sections[21]. 
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[_ Fig. 16 Computed near-surface streamlines for the initial and final HOPE configurations[20]. J 
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Fig. 17 Wing tip design of HOPE by KHI; initial and final chordwise pressure distributions 
over the tip fin and their spanwise sections[22]. 
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Fig. 18 USB wing planform[23]. 


Fig. 19 Sectional Cp distributions; initial and 
final Cp’s and their geometries[23]. 
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pressure distributions for initial geometry (upper surfee) 



Fig. 20 Wing-body combination design by FHI; initial and final chordwise Cp distributions and 
surface pressure contour plots[24]. 



Fig. 21 Final designed wing-body combination [24]. 
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